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We study a ghost-free model of massive vector curvaton proposed in the literature, where the quick 
decrease of the vector background expectation value is avoided by a suitable choice of kinetic and 
mass functions. The curvaton perturbations of this model have been so far computed assuming that 
these functions are external classical quantities, and it was found that some special time evolution 
of these functions leads to scale invariant and statistically isotropic perturbations of the vector 
curvaton. However, external functions should be understood as originating from the expectation 
value of some additional field. Since these functions need to present a non-trivial evolution during 
inflation, the field cannot be trivially integrated out, and, in particular, its perturbations need to be 
included in the computation. We do so in a minimal implementation of the mechanism, where the 
additional field is identified with the inflaton. We show that, except for a narrow window of model 
parameters, the interaction with this field generally causes the curvature perturbations to violate 
statistical isotropy beyond the observational limit. 
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I. INTRODUCTION 

In the past three decades, primordial inflation has become a dominant paradigm for the very early universe. It 
provides a simple, yet compelling mechanism to resolve the conceptual problems that the standard Big Bang cosmology 
confronts, namely the horizon, flatness, and monopole problems [TH3]- On the other hand, the observations of large- 
scale structure (LSS) and cosmic microwave background (CMB) exhibit small fluctuations of order 10~ 5 , which are 
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considered of primordial origin. The same theories that realize inflation can also predict such fluctuations with the 
observed features (near scale invariance, nearly Gaussian statistics, greater power in scalar modes than in tensor), 
which makes inflation a particularly favorable theory (see e.g. [I] for review). However, many realizations of inflation 
give degenerate predictions compatible with these observations, and therefore there is a large parameter space yet 
to explore for the underlying particle physics in the very early universe. Considerable amounts of study have been 
done to discriminate among otherwise degenerate models. One discriminator is non-Gaussianity, deviation from the 
Gaussian statistics, in the CMB. While non-Gaussianity is predicted to be undetectable in the simplest models of 
inflation [5-8 , a number of extended models have been studied which predict detectable non-Gaussian signature 
[5HT5]. The current bounds on non-Gaussianity are still loose [T5], but future missions, such as Planck satellite, are 
expected to improve measurement precisions to unprecedented level |17j . 

Observable non-Gaussianity, and other interesting signatures, can be obtained in the curvaton mechanism [18] . 
where the field responsible for the primordial perturbations is not the inflaton. While for simplicity it is most often 
assumed that the curvaton is a scalar field, models of vector curvatons have also been proposed [TM2"T] . Vector 
fields break the background isotropy through their vacuum expectation value (vev), and therefore their spectrum 
of perturbations may break statistical isotropy. Interestingly, a 30% level violation of statistical isotropy has been 
detected at 9a in the WMAP data [55J[23]. However, the direction of the asymmetry nearly coincides with the ecliptic 
one, which strengthens the case for a systematic origin of the effect. Fortunately, this can be checked by Planck, that 
can be sensitive to anisotropies at the few percent level |24j . The observed anisotropy, and the expected improvement 
from Planck, have motivated the study of cosmological models that can produce this effect, including those of vector 
curvaton. 

In the standard scenarios with an expansion by a (approximate) positive cosmological constant, it has been shown 
that the models of almost all Bianchi types evolve rapidly toward the de Sitter solution and any anisotropy present 
initially will be washed away (sometimes called cosmic no-hair conjecture) |25j . To obtain an observable anisotropy, 
the premises in [25] have to be violated. One way to accomplish this is to introduce a vector field with a non-vanishing 
vev during inflation. 

However, in the case of a familiar U(l) gauge field with a minimal kinetic term —F 2 /4 and no potential term, the 
field is conformally coupled to gravity. Thus its energy density quickly decays away and cannot source any appreciable 
anisotropy for sufficiently long time. By breaking the conformal invariance, the vector field can provide non-vanishing 
prolonged shear pressure, which supports a prolonged stage of anisotropic expansion and leads to violation of statistical 
isotropy. Several implementations of this idea have been suggested [5S1[37]; however, these models break the U(l) 
symmetry and inevitably introduce a longitudinal polarization which, for these models, happens to be a ghost [2"grf3"U] . 
To avoid this instability, a model was suggested in [3T], and further studied in [3"2"H3"4"] . where the U(l) gauge invariance 
is preserved and therefore no dangerous longitudinal mode is present. In this work, the fast decrease of the gauge 
field vev is avoided by introducing a scalar function that multiplies the vector kinetic term, L D — / (<p) F 2 /4, where 
ip is the inflaton field. 1 

A similar idea to that of [31] was used in [19l [20] , where the vector field plays the role of the curvaton, and where 
two time dependent functions were considered 2 

C = - f -^F^F» v - l -m 2 (t) A^A* (1) 

The kinetic coupling is typical of moduli or dilaton-like fields in string theory and supergravity frameworks, and the 
mass can be induced by a Higgs-like mechanism. Since the vector field is massive, the additional degree of freedom 
(longitudinal mode) is present, but, for /, m 2 > the model has no ghosts. It is shown in [THIEI] that the curvature 
perturbations generated in this model attain the scale-invariant statistically-isotropic power spectrum, if the varying 
functions have a specific time dependence, i.e. / cx a~ 4 and m cx a, where a is the scale factor, and if the vector 
mass is initially light and is heavy by the end of inflation (the physical vector mass m/y/J grows considerably during 
inflation). 

In a field theory, an external time dependence should be understood as the vev of a field. Given that the functions 
/ and m needs to present a nontrivial evolution during inflation, this field cannot be trivially integrated out, and, 
in particular, its perturbations cannot be disregarded. The field acts as "clock" ; therefore, the most minimal choice 



1 The kinetic term used in |31l is analogous to the one originally employed by Ratra |35| . The magnetogenesis application of I35| suffers 
of a strong coupling problem 1361 137 |; however it was shown in |37| that the model 1351 (without need of identifying the vector field with 
the electromagnetic one) produces non-Gauassianity of the primordial perturbations of nearly local type. 

2 More recently 1381 , this model was extended to includes the effects from a parity violating term, proportional FF, where F is the 
dual of the field-strength tensor. The non-Gaussian signatures due to particle production through the interaction ipFF, where tp is a 
pseudo-scalar inflaton, are studied extensively in |15| . 
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is to identify this field with the inflaton field. Ref. [35] already provided some specific examples where / and m are 
functions of the inflaton, and showed that, for a generic inflaton potential that satisfies the slow-roll conditions, the 
attractor solutions of the background evolution lead to small anisotropy in the expansion and to the required time 
dependence for the background values of / and m. However, ref. |39j does not compute cosmological perturbations in 
this set-up, but refers to the results of |19l I20j . in which the functions / and m are only external classical quantities. 

The non-minimal coupling through the kinetic and mass terms modulated by inflaton inevitably induces an in- 
teraction between the scalar and the vector perturbations. It is natural to think that this interaction modifies the 
evolution of perturbations in a non-trivial way. As we will show later in this paper, this interaction, which is present 
already in the linearized level, can actually be directionally biased due to the background anisotropy. Moreover, since 
the system of perturbations is now a coupled one, consistent quantization has to be done in the matrix form, the 



formulation first developed in |40j and summarized in Section III In this paper, we treat the full coupled system 
consistently, and show that the scalar-vector interaction produces direction-dependent effects in the power spectrum 
of curvature perturbations, in the framework of the vector curvaton scenario. We find that the near scale invariance 
of the spectrum is a generic feature of the model. The (near) statistical isotropy can also be achieved; however, it 
requires a specific choice of the coupling constant entering in the functions / and to, and appropriate initial conditions. 
This is in contrast with the approximated computation of [191 120) , in which the statistical isotropy appeared to be a 
generic result, independent of the functional form of / and to. 

We stress that, despite introducing an inflaton, we still want to work under the assumption of [19, 20 that this is a 
model of vector curvaton, so that the vector field should be responsible for the cosmological perturbations. This can 
for instance be achieved by assuming that, after inflation, the inflaton quickly decays into relativistic fields, and that 
the energy density of the decay products become completely negligible with respect to that of the curvaton (or its 
decay products). We do no follow the details of reheating here, but we simply compute the vector field perturbations 
until they freeze out, and we assume that they are the source of cosmological perturbations, according to the curvaton 
mechanism hypothesis. As also done in [191 120] . we still disregard (for technical reasons, as the computation is very 
involved) the metric perturbations in the analysis. These are the same working assumptions of the realizations of this 
mechanism suggested in |39j . In addition to |39j . we however consistently include the interaction with the inflaton 
perturbations induced by the two functions / and to that characterize this model. We show that this drastically 
changes the curvaton power spectrum in this model. 

The paper is organized as follows. In Section|llJ we analyze the background dynamics including the vev of the scalar 
and vector fields and the anisotropic expansion of the universe. The background attractor of this model is derived. 
In Section |III[ we discuss the evolution of the perturbations. We focus on the 2D scalar modes, which are the ones 
that contribute to the energy density and thus to the curvature perturbations. We quantize the system consistently 
for a coupled system and derive the equations of motion to evolve the system of perturbations in the matrix form. 
In Section |IV[ we numerically compute the power spectrum of the vector density fluctuations and relate it to that 
of curvature perturbations. Results are shown for some different values of parameters. In Section [Vj we discuss the 
results and conclude. Throughout the paper, natural units are used, h — c—1, and the Einstein notation is assumed 
for repeated indices. 



II. BACKGROUND DYNAMICS 



We consider a model with a scalar inflaton ip and a massive vector field whose mass and kinetic functions vary 
in time. This model is a realization of the mechanism of [191 120| , where cosmological perturbations were computed 
under the assumption that the vector field is a curvaton, and using the Lagrangian ([!]). The computation of [T51 |2"U] 
can be considered as an approximated computation of the perturbations, performed under the following non-trivial 
simplifying assumptions: (i) the background dynamics is in pure de Sitter, (ii) the background anisotropy induced 
from the vev of the vector field is negligible, (iii) the vector mass and kinetic functions are classical external functions 
of time, with no perturbations. The goal of this work is to improve over the results of |19L 120) by removing these 
assumptions. 



A. Model and Setup 



The model we consider is described by the action 



S = 



d A x- 



Ml 1 f(in) 1 „ 

' - -d^ d ^-V( l p)-^F fw F^--m 2 ^)A^ 



-R 



(2) 
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where F^ v — d^A v — d u A^ and M p is the reduced Planck mass. The universe expansion during inflation is driven by 
ip with a sufficiently flat potential V (p) , allowing slow-roll inflation (which we assume throughout this paper) . Here 
the vector mass and kinetic terms are function of p rather than being external functions of time. In [Tj|J HO] > it is 
claimed that the scale invariance and statistical isotropy of the power spectrum of primordial curvature perturbations 
can be achieved in the vector curvaton scenario, if / and m have time dependence / oc a~ 4 and m oc a, where a is 
the overall scale factor. With an appropriate choice of f (if) and m(p), these time dependences can be dynamically 
achieved (as shown in |39] , and as we review below for our particular model) . 

We can, without loss of generality, fix the background coordinate system such that the vector vev points in the x 
direction, i.e. (A^) = (0, A (t) , 0, 0). A background metric consistent with this choice is a Bianchi type I metric with 
the residual isotropy in the yz plane, given by 



ds 2 



-dt 2 + a 2 (t) dx 2 + b 2 (t) (dy 2 + dz 2 ) 



(3) 



where we parametrize the two scale factors by a (t) — e a ( t )- 2cr ( t ) and b (t) — Q a ( t )+ a '{' t ) . Hence a (t) is the number of 
c-folds averaged over all directions, and a (t) measures the difference in e-folds between the two directions. Note that 
a corresponds to the "overall" Hubble parameter, and a corresponds to the shear. 

There are four dynamical degrees of freedom in the background system of this model, cj> (t), A (t), a (t) and a (t), 
where 4> (t) is the homogeneous vev of ip. The corresponding equations of motion are, respectively, 



3qi 



V' = fe 



A 



a - 



46- 



r 
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3& 2 = 



M 2 



/ e -2a+4er 



6M p 2 



/ 



-2q+4ct 



3M p 2 



A 2 



-A- 



with the Einstein constraint equation 

3d 2 - 3cr 2 



1 

M 2 



<t> 2 + V ) + 



-2a+i<y 



\2 



A 



-A 2 



f 



-A* 



(4) 
(5) 
(6) 
(7) 

(8) 



Here, and throughout the paper, dot denotes the derivative with respect to t, and prime with respect to p. Notice 
that Eq. ([T]) implies that the evolution of a is directly supported only by the pressure of vector field. Also notice that 
the system is invariant under the unphysical rescalings a. — > a + ao, a — > ct + cto, A — > Ae Q °~ 2 °' , where ao and cto are 
arbitrary constants, and thus the normalization of a and a does not affect any physical quantity. As a consistency 
check, if we set A = and a = 0, we would recover the standard single-scalar-ficld inflation. It is also useful to define 
the background energy densities for the scalar and vector fields, given by 



PA 



1 
2 
f 



V 



-2o+4a 



.4 



i2 



m 

T 



A 2 



(9) 
(10) 



The inflaton energy density is the standard one, and in the slow-roll regime, V (s> \4> 2 ) drives the quasi exponential 
expansion. Due to the non-zero vector mass, the vector energy density is also split into the kinetic and potential 
parts. Interestingly, in the pure de Sitter (& = and a is constant), pa is exactly constant during inflation, if / oc a~ 4 
and m oc a. As a consequence, in the quasi de Sitter background with small anisotropy, which is the case we consider, 
Pa stays nearly constant during inflation until the desired time dependence of / and m starts to be violated near 
the end of inflation. Also, it is worth noting that A is not a physical quantity; it always appears in the combination 
fe~ a+2a A, which is the physical one, and the physical mass of the vector field is 



M 



(ii) 



In the following, we occasionally use M instead of m, when the meaning is more transparent. We should note that 
M scales as oc a 3 during inflation and that we are interested in the parameter space where M <^ a at early stages 
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of inflation and M ^> a by the end of inflation, which attains the scale invariant and statistically isotropic power 
spectrum in the curvaton scenario for the case of isotropic de Sitter background and unperturbed /, to |19| 120] . 

The desired time dependence of the mass and kinetic functions, / cx a~ 4 and mtxo, have to be achieved dynamically. 
For a concrete realization of this, we consider the simplest chaotic potential of inflaton 

V^)= 1 -m 2 y. (12) 

With this form of potential, in order to obtain the desired time dependence, / and m take the form 

2 2 

f(<p) = e^, m(^)=m e" d % (13) 

where c and too are positive constants. Notice that after inflation, <j> (t) starts oscillating and approaches 0; accordingly, 
/ approaches to unity, and to to a constant value tuq. We will now derive the attractor of the background system 
and show that it leads to / cx a~ 4 and to cx a during inflation. 



B. Attractor Solutions 



We are interested in finding the attractor in the slow-roll inflationary regime. There are four dynamical background 
degrees of freedom in this system, but none of them has an exact analytic solution. However, in the pure de Sitter 
with negligible anisotropy, (pSJ) has an analytic solution for A (t) in the case / cx a~ 4 and to cx a. Motivated by this 
fact, we parametrize A (t) and A (t) in terms of C (t) and D (t) in the following way: 



A(t) = C (t) cos J dt'M [0 (t')] + D (t) sin J dt'M [cf> {t')\ 

-C (t) sin dt'M [(/> (t')] + D (t) cos jf dt'M [0 ft')] 



A(t) = M[cj>{t)} 



(14) 
(15) 



where t- m is some early time during inflation. Such parametrization can always be done. In pure de Sitter, C and 
D are integration constants; here, they arc promoted to be time dependent. In this parametrization, we can solve 
two first-order differential equations for C and D, instead of one second-order equation for A. The former is in fact 
numerically favorable, since C and D vary little during inflation and do not have a large hierarchy between themselves. 
Also, we assume the initial equipartition for the vector energy density: namely, the vector kinetic energy is equal to 
the vector potential at initial stages of inflation (this is assumed in [T!5], and we verified that it is indeed a necessary 
condition for the statistical isotropy in the power spectrum in the case of unperturbed / and to). Since M <C a at 
early times, we have A ~ C and A ~ MD, and thus the energy equipartition implies 



(16) 



To investigate the background attractor solution, we assume |C| ~ \D\ and that C and D are constants (the precise 
solution of C and D is instead used in the cumputation of the perturbations). 3 

We focus on the background equations for 0, a and the constraint, namely (Ej), ^ and Q, respectively. Here we 
consider the early time regime, namely the vector mass is light M -C a. From (Wj), ignoring cf> and a terms (slow roll 
and small anisotropy), we have 



3 a i 



Now we use a as the time variable with 4i 



l 2 (C 2 + 2D 2 ) -3a- 



2 Ml 



(17) 



will be verified), we have a 2 



3 Ml 



V. Defining 



4il/; 



ot~r~- Since the right-hand side of \d>\ is dominated by the term V (which 



da 



-2M„ 



■/?) 



a, we write the above equation in terms of . 
(C 2 + 2D 2 ) e"^%. 



(18) 



3 In fact, D changes some, depending on the value of c and the value of a (which in turn depends on the ratio between the vector kinetic 
and potential energies), while C stays really constant until near the end of inflation. 
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FIG. 1: Time dependence of /. The z-axis is the overall e-folds a, and y-axis is the quantity f/ fa, which would be equal to 
—4 if / oc a -4 exactly. This desired time dependence is well retrieved until the vector field becomes heavy (a ~ —2.2). At 
M ~ d, the time dependence is slightly violated; before it comes back to —4, inflation ends (a = 0), and / approaches 1. The 
parameters here are taken to be c = 1.5 and mo = 1000 m v . By construction, m obeys oc a when / oc a -4 . 



The homogeneous solution for this equation implies that the time variation of 
left-hand side of (18 1. Then we find 



c 4> A 
2 M 2 



c-1 IMl 



in. 



c 2 C 2 + 2D 2 ml 



is small, so we can neglect the 



(19) 



Since the right-hand side of ( 19 ) is constant, we obtain a 



- 4 M 2 , up to some constant. This means that the chosen 
forms of / and m indeed obey the desired time dependence, / oc a~ A and m oc a, in this attractor. This is verified 



numer ically in FIG. [IJ Combining this 



s with 



l 

6 Ml 



m 2 cf) 2 gives 



3 a < 



m%4> 



m v M p 



(20) 



This coincides with the attractor found in [3T]. It is worth noting that for this attractor solution to be valid, we need 
to have c > 1, as can easily be seen in ( fl9| ). From here on, we assume this range of values of c. 

To consider the anisotropy, we ignore the a term in ^ and take the phase in the parametrization of A and A to 
be zero. Since we assume initial equipartition for the vector energy density, it is convenient to define 



PL 



(21) 



i.e. the ratio of the vector potential to the vector kinetic energy density. We are interested in the case where this ratio 
stays close to 1. At initial times, this ratio reduces to r™ tlal 
we find the attractor lead the anisotropy to 



D 2 



. Combining these with ( 19 1 and a 



6 Ml 



a 

a 



4 c 
3 ~< 



lr A -l M 2 



(22) 



Notice that a vanishes if ta = 1- This is clear both intuitively and from ^7j. If = p^ 111 , the right-hand side of 
becomes 0, and then a decays away to quickly. The only matter content in this model that causes the shear pressure 
is the vector field. When the vector kinetic and potential terms arc equal, the pressure vanishes, like non-relativistic 
dust, and does not enhance anisotropy. In the actual evolution, ta does not stay exactly equal to 1, so anisotropy 
is produced to a small amount, until the vector field becomes heavy. Once the vector mass exceeds the background 



7 
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FIG. 2: Evolution of the anisotropy, |<r|/d. The numerical solution (blue solid curve) is compared to the attractor solution 
(22 l (red dashed curve), where the numerical solutions are used for ta and <j>. Here the overall e-folding (a) is normalized to 
when inflation ends. The vector field becomes heavy (M = d) at a ~ —2.2. As discussed in the main text, the attractor 
solution matches well with the numerical until M w a; after this moment, the anisotropy decays quickly, but (22 1 does not 
replicate that. Here the parameters are taken to be c = 1.5, mo = 1000 m v and r^ 1 * 13,1 = 1. For other values of parameters, 
the attractor is valid as long as c is not too close to 1. 



expansion rate (M > a), the pressure vanishes and so does the anisotropy. FIG. [2] shows that (22 1 is a very good 
approximation, until the vector field becomes heavy (M ss a). 

We can also determine from the attractor the ratio between the energy densities of the scalar and the vector. 
Assuming slow roll and small anisotropy, we find 



PA 



IrA + lMl 

r A + 2 <j> 2 



(23) 



We note that in the limit of ta — > (massless), this recovers the attractor solution found in [3TJ [33] and one of the 
attractor solutions in [39]; on the other hand, in the limit of ta S> 1, it recovers another attractor in [39]. We stress 



that our early-time attractor expressions ( 20 1 , ( 22 ) and ( 23 ) are valid for arbitrary values of Ta 



Our starting assumptions are now verified: anisotropy is small {a <C a), and the vector energy density is subdomi- 
nant (pa <C p v ) , in the slow- roll regime with the potential V of the form ( fl2| ) . The vector kinetic and mass functions 
of the form ( 13 ) are dynamically led to their desired time dependence by the attractor. The existence of the attractor 
requires c > 1, and we only consider this range of c. Also, we assume the initial equipartition for the vector density, 
which translates to \C (ii n )| = \D (ti n )| as an initial condition. As long as anisotropy is small with quasi de Sitter 
expansion and / and m obey the time dependence, C and D are nearly constant, and the equipartition is maintained. 
In fact, if c is too close to 1, the attractor is no longer a good approximation; then the time dependence of / and m 
starts differing from the desired one, and the initial equipartition deviates as the system evolves in time. 



III. PERTURBATIONS 



We now consider the perturbations of the model (|2j). To avoid an extreme computational complication, we neglect 
metric perturbations in this paper. While a further study including them would be interesting, we will show that the 
inclusion of background anisotropy and the perturbations of ip drastically changes the power spectrum of curvature 
perturbations, as compared to the case of isotropic, pure de Sitter background with no Sip, studied in |19| . Here we 
study the effects in the linearized level of perturbation theory. 



A. Basics 



First note that since the vector mass is put by hand, there is no U(l) gauge freedom (so we never call it "gauge" 
field). Consequently the vector perturbations contain three independent degrees of freedom, two transverse and one 



longitudinal, in terms of momentum directions. In some of massive vector models, the longitudinal modes may suffer 
from ghost instabilities [28Tl30| : however, this particular model we consider here is ghost-free [20] . 

Since we are neglecting the metric perturbations, the metric still takes the form ([3|. We define the field perturbations 
as the field values with their background subtracted, 



and their Fourier transformations through 

d 3 k 



Sip (t,x) 



(27T) 



3/2 



SA^ (t, x) = A„ (t, x) - (Ap) 



e***M*,fc , tiA IM {t,x) = 



d 3 k 



(27T) 



3/2 



e lk ' x SA„ t, k 



(24) 



(25) 



as usual. To simplify the computation and to make geometrical interpretation transparent, we orient the coordinate 
system such that the z-component of the comoving momentum vanishes |41j . 



k = [k L ,k T ,0) 



(26) 



where L and T denote "longitudinal" and "transverse" in terms of the background vector vev. This choice can be 
made without loss of generality, due to the residual 2D symmetry in yz plane after orienting the vector vev along 
the IE-axis (as done in Section [IT]) . In this coordinate system, it can easily be seen that the 2nd-order action is split 
into two decoupled sectors, 2D scalar ({Sip, 6Aq, 5A x , tiA y }) and 2D vector ({SA Z }). The 2D vector mode 6A Z does 
not contribute to the curvature perturbations, which are intrinsically a scalar quantity (this is more explicitly shown 
in Section IV), and thus we disregard this sector in the following discussions. We denote the physical momentum 



p = (j>l,Pti 0): with its components and norm 
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(27) 



(Recall a = e a - 2a and b = e Q+<T .) 

Among the 2D scalar modes, SAq is non-dynamical and can be integrated away. We see that SAq is related to other 
2D scalar modes by 



6An = 



M 2 
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8A X 



^SA 
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f a 



(28) 



We now have three independent modes in 2D scalar and one in 2D vector, as expected. One of the two vector 
transverse mode is the 2D vector mode SA Z ; the other transverse and the longitudinal modes are linear combinations 



of 5A X and 5 Ay. The 2D scalar sector of 2nd-order action is now 
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(29) 



We need to solve this system to find the power spectrum of curvature perturbations. However, it is not possible to solve 
such a complicated system exactly; moreover, we look for the effects from non-zero a and from the coupling between 
SAfj, and dip, and thus approximated methods would be quite delicate. Thus we employ numerical computation here 
to obtain the final results. To this end, we need to appropriately quantize the system and find the initial conditions. 
This is a coupled system, and so the quantization and the evolution have to be done in the matrix form, as described 
in detail in |33l 140] and summarized in the following subsections. 
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B. Quantization 



We first rotate the system such that the kinetic mixing term present in the first line of (29) vanish. This can be 
done through 



Si T&ijYj 



where Si = (Sip,5A x ,SA y ) i and a rotation matrix 



Hij = e" 3Q / 2 



/I 












■■ C T 



— 



'■r r=, ,„ c l 



(30) 



(31) 



with arbitrary time dependent functions ct and cj,. By choosing them as 

p 2 + M 2 



C-rp 



2/ 1+plPt/V(pI+M 2 )(p 2 t + M 2 ) 



cl 



p 2 + M 



2f l-p L p T /y/(p 2 L + M 2 )(p 2 T + M 2 ) 



(32) 



then the kinetic matrix becomes the identity. Subtracting appropriate total derivatives, the action (29) takes the form 

1 



o(2) _ 
°2dS — 2 



dt d A k 



yt y + yt x y - yt x y - yt n 2 y 



(33) 



where all the components of the matrices X and Ct 2 are real. These matrices have further properties X T = —X and 
Cl 2T = fl 2 , which is always possible by subtracting total derivatives; moreover, X and Cl 2 are invariant under the 



parity k — > —k. Also note that all the field components entering in the action satisfies Y, 



l(-k)=Y(k) ( 



so that y 



are real in the coordinate space) . The action ( 33 1 is formally of the same type as in [33 , so we follow the procedure 
outlined there. 

To remove the X terms, define 



tp = RY 



(34) 



where R is a real, orthogonal matrix, so that Y^ Y = ip^tp, and satisfies the relation R = RX. The action now 
becomes 



S 



1 



dtd 3 k 



^tp-^n 2 tp 



(35) 



where Cl 2 = R (Cl 2 — X 2 ) R T . Note that we can always choose R to have the property R (—k) — R (k) since 
X satisfies the same property. Hence, f2 2 \—k) = fl 2 (fc), and tp^ \—k) = tp (k). Thanks to this property it is 



immediately realized that tp is an array of real fields in the coordinate space and that the conjugate momentum is 



defined as iii = tp\. The Hamiltonian is then H = | J d 3 k it^tt + ip^Srtp 
define 



1p = C tp , 7T = C 7T 

where C is an orthogonal matrix (C T C = 1) that diagonalize £1 2 , i.e 

C T Q 2 C = diag (w 2 , uj 2 , uj 3 

giving H = \ J d 3 k 7T+7T + ■0 t w 2 ^ . Note that, since f2 is a Hermitian matrix and O 



To diagonalize the Hamiltonian, we further 

(36) 



• 2 ^ - u 2 



(37) 



-k) = 



(jtj . Since C is unchanged under the parity k — > — k, we then define 



(I) 



, uif are all real 



P* «t 

1 lm u m 



tp (fc) = Cy (fc) [fy, (fc) P lm a m (fc) + h* x (-t) 

(k)a m (k)+h* l (-k)PLal(-k) 



7T ( k 



k\Pi, 



(38) 
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where a, (jtj and aj ^fc^ are arrays of annihilation and creation operators, respectively, satisfying the commutation 

relation cij ' a I (^ 2 ) = — ■ ^ ne ma trix Py is a constant Hcrmitian matrix; this phase freedom is 

present in quantizing a coupled system, similarly to arbitrary phases appearing in a free-field theory, and this phase 
does not affect any physical result (in fact, any consistent quantization should be constructed allowing such freedom). 
In the coordinate space, tp and tt are real, and to have ip and tt satisfy the equal-time commutation relation, we require 



h k)hUk) -h* -k) h T (- 



7 (-i) 



= i Si 



We decompose h and ft in a similar manner to Bogolyubov transformation: 

1 „^ r —iuJ 



h 



'2lo 



■0) 



'2uj 



(a-p). 



(39) 



(40) 



Here h, h, ui, a and P are all matrices, but since ui 2 is a diagonal matrix, l/v2u) and w/v2w are well defined matrices. 
With these "Bogolyubov coefficients" a and /?, the condition (39) is translated to 



a (jfc) a f (it) - p* (-k) p T (-lej 



a (k) ft (fc) - 0* (-k) a T (-k) 



0. 



(41) 



This is not the unique choice to satisfy ( 39 1 , but it is a sufficient condition and the following derivation is consistent 

(42) 



with this choice. The Hamiltonian now becomes, after normal-ordering, 



where 6, and b] are defined as 



b[k 

&t (-k 



:H:= I drkwib] 



a ( k 



k) b„ 








P I k) a* [ —k 



a (it) 
a) (-% 



(43) 

The Hamiltonian is thus fully 



Using (41), it can be shown b and satisfies bi (ki) , &j (^2) = % 5^ (jti — k 
diagonalized and the system is consistently quantized. 

We are now ready to find the equations of motion for the system in terms of h\ or of {a, P}. As we have 

ipi (jtj = -Ki (jtj and -hi (jtj = —Q 2 (jtj ip (jtj (from (35)), the equations of motion for h and h are 



h = h — r h , h = — r h — ui 2 h 



where T = C C. With a and /3, these equations are expressed as 

to 



a 



-i to a H P — la — J P 

2u> 



P = iuj P + — a-ip- Ja 
2ui 



where 



fuT- 



J 



TJZj 



(44) 
(45) 
(46) 



At asymptotically early times, when modes are deeply inside the horizon, flfj ~ p 2 5ij and Xij <C p, and consequently 

<jji ~ p and r, J, J, — -C oj. Thus we obtain the adiabatic initial condition a- ln — e~ l $ dt Win l and p- m = 0, or 
equivalently 

1 _j r* At',.,. 1 



-i /* dt'u 
V 2 w in 

ln — m 



3 -t/' dt'u tB j 



(47) 



This is the solution to ( 44 ) at early times and satisfies the quantization condition ( |39[ ) . It corresponds to the initial 
adiabatic vacuum, if we choose the Bunch-Davis initial vacuum which is annihilated by the operator ai appearing in 
p8| . Therefore, we have formally quantized the coupled system of the 2D scalar sector of the model, and we can 
determine the initial conditions through ( 47 1 to evolve the system. 
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C. Early-Time Evolution 



We now need to evolve the system given by the action (29), with the initial conditions determined by (47). As is 



clear from the quantization procedure, this must be done in the matrix form. It would be ideal to rotate the system 
using the rotation matrix (31 1 with the coefficients c t,l in (32) to make the kinetic matrix the identity, and to evolve 
hij and hij according to their equations of motion (44). However, we have found it computationally challenging to 



use (132). To circumvent this issue, we still rotate the system by (31) but choose ct,l such that the kinetic matrix 
becomes very close to the identity only at early times M<Cp. For this purpose, we take 



c T 



P 



2v7 



CL 



PL PT 

A/v 7 /' 



The action ( 29 1 now becomes 



o(2) _ £ 
°2dS — 2 



dtd k 



r f ty + y 1 xy + y 1 x t y - r f n 2 Y 



(48) 



(49) 



Here, each component of Yi is related to that of Yi in the previous subsection only by rescaling due to a different 
choice of the coefficients. The matrix T is diagonal, and at early times, we have T ear i y ~ 1 as we required; moreover, 
^cariy — P 2 1 an d ^cariy Pi necessary for adiabatic initial conditions. The complete expressions for the matrices T, 
X and O 2 are written in Appendix[A] From the expression of T33, it is manifest that I3 corresponds to the longitudinal 
mode, with the free choice of sign (plPt) > 0. 



We evolve the system of Yj. By varying the action (49), the equations of motion are obtained as 



where 



Y 



■A-ij Yj ~h $ij Yj 



(50) 



A 



T- 1 (f + X- X T 



B = T~ l fi 2 + X 



(51) 



Since T is diagonal, T 1 is trivially found. We now need to decompose Yi into arrays of creation and annihilation 
operators with matrix coefficients, consistent with the quantization procedure, as in (38 1. That is, 



Y 



(t, kj = y i3 (t, kj a, (kj + y*j (t, -kj a] (-kj 



(52) 



where the creation and annihilation operators here are the same as those that enter in ( 38 ) . The matrix coefficients 
are related to hij by 



y, 



(r-l/2 R T ( 



Ch) P, 



Ij 



y, 



T- 1/2 R T Ch + d t (r- 1/2 R T cj h 



P, 



ij 



(53) 



where the mat rices R, C, and P are introduced in the previous subsection. By writing in terms of y, the equations 
of motion (50) consist of two parts, one proportional to a 3 (kj and the other to (—kj , an d these two terms must 
vanish simultaneously. They are in fact not independent but equivalent terms, due to the fact that A and B are real 
and invariant under k — > — k. Therefore the equation of motion we need to evolve in the matrix form is 



y., + An yij + Bu y>, = 0. 



(54) 



To determine the initial conditions, we employ the adiabatic initial condition (47) and relate it to y through (p3 



In doing so, we take advantage of the phase freedom by fixing the arbitrary constant Hermitian matrix P to eliminate 
the initial phase. (This procedure is done in [31].) Namely, we choose P — C T R at initial time; since both of R 
and C are orthogonal, and since both h and h are proportional to the identity, we can simplify (53) to y = T^ 1 / 2 h 
and y ~ T ' h + na T~ 1 / 2 h, where n is some number of O (1). Recall that T oar i y ~ 1 with subdominant terms 



4 Here, X is not an antisymmetric matrix, i.e. X T 7^ —X. Although anytisymmetrization is always possible for X, we did not do it 
here, only for simpler expressions of the matrices, which are explicitly written in Appendix|A] This does not change the evolution, since 
antisymmetrization is done by adding total derivatives in the action and thus does not affect the equations of motion. 
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completely negligible (suppressed by ~ M 2 /p 2 ), and d t (T x / 2 R t C) ~ aT 1 / 2 R T C. The second term in y is 
suppressed by a/p and thus is negligible. Therefore, we give y initial conditions 



1 

\/2pin 




in — /rf-^ — ) j, in — M/ o ^y- (^5) 



We evolve the system in terms of y at the early stage of inflation, i.e. until a few e-folds before M = p. 5 After 
this moment, we switch back to the original variables Si = {5Lp,SA x ,5Ay), for the ease of computation of the power 
spectrum. In the next subsection, we briefly formulate this late-time evolution. 



D. Late-Time Evolution 



For the late time of evolution (after M ~ p), we rotate the system back to the original variables Sip, SA X , SA y . We 
do this as it is numerically efficient and the computation of power spectrum is much more compact. We could not do 
this for the early time, due to a huge hierarchy between the transverse (corresponding to Y2) and longitudinal (to Y3) 
modes of the vector field. This hierarchy is originated from the fact that M -C p at the beginning (to have a feeling, 
if we take tuq = 1000 m^, then Mjp <~ 10~ 76 at horizon crossing). The fields SA X and 5A y are linear combinations of 
these modes, and no sensible numerical precision can take such tiny initial values into account. This is why we had 
to evolve y for early times. After M > p, however, the transverse and longitudinal modes evolve similarly, and thus 
such an issue no longer arises. 

The action for the late time takes the form 



;(&TS + tf XS + S^ X T S + S^n 2 s"j (56) 



where Si = (Sip,SA x ,SA y ) and T is not a diagonal matrix. This is nothing but a compact expression of (29), and the 
matrices T, X and ft 2 can easily be identified. As for the early time, the late-time variables must be evolved in the 
matrix form. We decompose Si by 

Si (t, tj = Ay (t, fc) aj (fc) + A*j (t, -fc) a] (57) 



where the creation and annihilation operators are as in (38) and (52). At the transition, we connect the late-time 
matrices to the early-time ones by 

\j - Tin yij , Aij = TZa y tj + lla y h . (58) 



The equations of motion for A^- can be found in the manner similar to the early-time case; varying the action (56), 
identifying the terms proportional to aj (jtj and ^— kj , we obtain 

+ A a Ay + B a Ay = (59) 

where 

A = f- 1 (f + X- X T ^j , B = T- 1 ((l 2 +xy (60) 

Using the numerical solution for Ay, we compute the power spectrum of primordial curvature perturbation in the 
vector curvaton scenario. In the next section, we formulate the procedure and show the numerical results of spectrum. 



In numerical computation of the early stage of evolution, we needed to Taylor-expand each component of all the matrices up to sufficient 
order. Without doing this, the components that are relevant for the longitudinal mode of the vector field encounter cancellations in 
the leading order in the expansion. The sub-leading terms are, however, proportional to M 2 /p 2 , which is too tiny at early times to be 
within numerical precision. Such cancellation can be seen in the matrix components T33, X$i and in Appendix lAl 



13 



IV. OBSERVABLES 



A. Vector Curvaton Scenario 



The mechanism we consider in this paper to produce the curvature perturbations is the curvaton model. It is the 
mechanism to convert isocurvature (entropy) perturbations into curvature perturbations [TSJ [55] (see also [23J03]). In 
this model, the background expansion is driven by the inflaton (</>, in our case), but the primordial curvature pertur- 
bation is generated by some other field "curvaton" (A^), while the direct contribution from the inflaton perturbation 
on the final observed cosmological perturbations is negligible, unlike single-field inflation. However, even if the model 
is arranged such that at late times the contribution of the inflaton perturbations to the cosmological perturbations 
is much smaller than that of the vector curvaton, the inflaton perturbations can still drastically affect the vector 
perturbations due to the direct ip — A^ coupling that necessarily originates from the two functions / (if) and m (if). 

The physical picture we have in mind is the following: the primordial perturbations of both fields if and A^ are 
originated from quantum fluctuations in the initial adiabatic vacuum. At some time during inflation, the vector 
curvaton becomes heavy and starts oscillating; however, due to the coupling with inflaton that breaks the conformal 
invariance of the vector field, the vector energy density does not decay away even after the oscillation starts but 
stays nearly constant until the coupling is terminated at the end of inflation. Some time after inflation ends, the 
inflaton decays into radiation (the detailed mechanism of (p)reheating is beyond the scope of this paper, and we 
simply assume an instantaneous reheating not long after inflation). The curvaton still oscillates, but it now behaves 
as non-relativistic dust, since its mass is larger than the Hubble parameter and its coupling to inflaton is no longer 
present. Since the radiation energy density decrease as oc a~ 4 and the oscillating curvaton density as a~ 3 , the 
cosmological perturbations induced by those in the radiation density becomes negligible. Due to the coexistence 
of radiation and dust, the pressure is non-adiabatic; as the curvaton oscillation continues for many Hubble times, 
this non-adiabatic pressure perturbation converts the curvaton perturbation into the curvature perturbation. Finally 
(and before neutrino decoupling |18j). the curvaton decays and the curvature perturbation remains constant until the 
horizon re-enetry. We assume for simplicity that the curvaton decay is also instantaneous. 

We consider the curvature perturbation after the universe becomes isotropic (which is true after the vector field 
becomes heavy), and work in the spatially flat slicing. In this slicing, the curvature perturbation defined on the 
uniform-density hypersurfaces is given by [181 145] 

c = -^^zrr— (61) 

p 4 + 3 r p A 

where p = p r + pA is the total background energy density, r = PaI Px is the ratio of the background vector energy 
density to the radiation density, and bp a corresponds to the vector energy density perturbation. In the approximate 



equality in (61), we assume that the radiation perturbation has already become negligible at the time of curvaton 

1 SpA . 
3 pa ' 

_ In on -f- Vi err nocci -fVid I ocoAn Vi c\ra it: f ha+ fniD l'mnArTOTit niicmfifir uro naoH + r\ rr\ r» 1 1 ■ 

4 PA 



decay and that p r oc a 4 and pa oc a 3 . If the curvaton dominates the energy density before its decay, then £ ~ - ^ 
in the opposite case, £ ~ 5 177"- ^ n either case, the lesson here is that the important quantity we need to focus on is 
the ratio 5 pa/ Pa- 

We can identify the energy density as — Tq, where T£ is the energy-momentum tensor of the model, and its 
perturbation is found by perturbing the field contents. In the limit Sg^ — 0, we then find 



(> 



SpA = 

a 



f 



' + M bp 2 + M 2 y 



P 



(62) 



in the Fourier space. In finding (62), SA is already substituted by (28). Notice that (i) the unphysical normalization 
of scale factors and A^ does not affect Sp, which is a physical quantity (see the brief discussion after (|8|), and (ii) the 
2D vector mode SA Z does not contribute at linear level to 5p, which is intrinsically a scalar quantity. This justifies 
disregarding the 2D vector mode. 



The 2D scalar sector is a coupled system and the quantization is done with respect to the initial vacuum, so, in principle, also the 
cross-correlation {5pA$P<p) contributes to the value of (8p5p). We verified numerically that the cross correlation is comparable to 
(&PA a) a t the end of inflation. Since our focus is the study of a curvaton mechanism, we want to assume that &p v plays no role in the 
post-inflationary perturbation. This can be trivially achieved if, after inflation, the inflaton decays into a relativistic species X before 
the massive vector decays. Then there is a stage in which px oc a -4 , while pA oc a -3 , and the cross-correlation decreases a factor of a 
faster than (&Pa 5p<p)- 



14 



By using the decomposition ( 57 ) , we write 



bp A (t, k) = 5p A ,i (t, kj a,i (k) + 5p* Ai ( 



t, -k 



SpA,i 



I 



p2T + M2 AA 2l - a T ^,AA 3l 



4H 



_ p 2 +M 2 « b P 2 +M 2 
The 2-point correlation function of curvature perturbation is (£ 2 ) cx (Sp A ) / 1 Pa-, so we compute 

1 ,„ , _ „ . ^, /" d 3 fc 



7 0?, y) 



(5pA (*, x] bp a (t, y)) 



e ik-(x-y)-p 



where the expression of is shown in (10), and 



(63) 
(64) 

(65) 
(66) 



When wc chose the momentum orientation in (26), we set the coordinate system so that k 3 vanish. To be consistent 



with this choice, we set x — y = (tl, Tt, 0). Since the function T is unchanged under the parity k — > — fc, (65 1 becomes 

dk A 



7 = 



d£ cos(fc£r L ) J (k y/l - e r T ) 



(67) 



where k = \Jk\ + fc 2 -,, £ = fcz,/fc, and J v is the Bessel function of the first kind. We have introduced the power 
spectrum 

P = Jt 
7 " 2^ 2 



(68) 



Notice that, in the case of statistical isotropy, T depends only on k, and we recover the standard expression 7 

3st is 

p 2 T + M 2 



I IT ^"fcf" 1 ^7 ■ Finally the quantity of interest is 



„ _ fc 3 P ^ 

T ~~ 9^2 „4„2 



2tt 2 a^p\ 



M 2 



A A 



2 1 



a p L Pt 
b p 2 + M 2 



A A 3 j + M A A 2 * 



and in the late-time limit, the middle term in (69) is completely negligible and p <C M, so we have 

2 



■ 7 



2fc 3 

7T 2 



E 



A 2 + M 2 A 2 



(69) 



(70) 



P 7 . 



where we used (10). This is the quantity we compute. In the following subsection, we show the numerical results for 



Due to the non- vanishing anisotropy in the background expansion (until the vector field becomes heavy) and due to 
the scalar-vector interaction, it is natural to consider the possibility that a directional dependence of P 1 arises, while 
it is absent in the simplest single-field inflation. To quantify this statistical anisotropy in the spectrum, we employ 
the ACW parametrization 



P (k) = Piso (fc) [1 + g* e] 



(71) 



where £ is the cosine of the angle between the mode fc and the background privileged direction, coinciding with the £ 
abo ve. 7 This parametrization shares the same concept as the multipole expansion for small higher-order moments. In 
( 71 1, parity (fc — > — fc, £ — > — £) symmetry is assumed, which is valid in the model of our interest. Therefore the leading 
effect of anisotropy is proportional to £ 2 , corresponding to a quadrupole. Let us remind that this expansion is valid for 
small and for |^| ~ 1, the contributions from higher-order multipoles become important. The detected statistical 
anisotropy that we have mentioned in the Introduction is = 0.29 ± 0.031 [22 [23], although, as we mentioned, 
this effect is likely systematic. As a reference, we work under the assumption that the primordial perturbations are 
statistically isotropic with \g t \ < O (10 _1 ). Planck will probe the anisotropy up to \g* \ =0 (10 -2 ), as shown in |24j . 



In the case of bispectrum, the parametrization of statistical anisotropy in a vector field model is developed in I46| . 
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FIG. 3: Power spectrum P 7 denned in (70 1, in the units oimt/M^. The left panel shows the result when mo = 50 m v , and the 
right when mo = 1000 m^. In each plot, both the complete results (solid curves) and the cases where the coupling is artificially 
neglected (dashed curves) are shown. For each case, a few values of £ (corresponding to different orientations of p) are taken 
to illustrate the directional dependence. Here c = 1.5 is arbitrarily chosen to clearly exhibit the difference. 



B. Numerical Results 



In this subsection, we show the result of the power spectrum of primordial curvature perturbation (up to overall 
factor 1/3 or r/4) computed by (69 1 (or its late-time value (70)). As we discussed in Section [TlJ the initial equipartition 
of the vector energy density is required for statistical isotropy in the de Sitter limit with 5p = [15], i.e. fy n = p\, or 
|C (^in)| = \D (tin) |- We take this condition for all the numerical computation. For the same reason, the vector mass 
needs to be M -C a. initially and M 3> a by the end of inflation, and thus we choose the values of too accordingly. 
Also, to ensure the existence of the attractor ( [20[ [22] [23] ), c > 1 is required. We concentrate on the parameter space 
that satisfies all these conditions. In numerical calculation, we take £ = Pl/p evaluated at the end of inflation. This 
coincides with the definition in the previous subsection, since the vector field is already heavy at the end of inflation 
and the expansion is isotropic from then on. 

Our main interest is the effects due to the non-minimal couplings between the vector and scalar fields introduced 
in the action ([2]). These couplings are controlled by two parameters c and m , so in order to see such effects, we 
observe how the spectrum changes for different values of these parameters. FIGs. [3] and [4] show the computed power 
spectra for a few different values of £. 8 These choices of £ are solely of illustrative purpose and have no particular 
physical significance. In this model, the coupling functions / and m smoothly approach to constants after inflation, 
but this does not occur instantaneously. As a result, it takes some time for the system to completely decouple and 
for the curvature perturbation defined in (61) to become adiabatic direction and freeze out. We therefore evaluate 
the spectrum a few oscillations (of inflaton) after the end of inflation, when the variation of the values of P 1 is at 1% 
level. 

FIG. [3] shows two cases that compare different values of vector mass: too = 50to v on the left panel and mo = 
1000 m v on the right. In each panel, two plots are compared, one (solid curves) with the consistent evolution of 
the coupled system, and the other (dashed curves) with the vector-scalar coupling neglected for the perturbations 
(equivalent to setting Sp — 0). Let us emphasize that in the latter, setting Sip = is completely artificial and is 
inconsistent within the model. We show this unrealistic case only to illustrate how significant effects the coupling 
produces, and one can observe that the impact is indeed substantial. If the coupling is manually turned off (dashed 
curves), the observed or smaller statistical anisotropy can be attained, consistent to the claim in [19] [20] 9 ; however, 
if the coupling is present, then the statistical isotropy is clearly violated (g* ~ —0.8), as can be seen in the figure. 



The horizontal axis of FIG. [3] and [4] is po/ag ~ k/ (aoifo), where denotes the end of inflation. Ideally, we would compute the spectrum 
of the modes that correspond up to the largest observable scale, or the modes that exit the horizon 60 e-folds before the end of inflation, 
but due to numerical difficulties, we show the modes pq/cxq = 10 -10 to po/ao = 10~ 3 and expect no significant change in the spectral 
behavior for smaller momenta. 

In fact, this is true only for the coupling constant c > 1.25. If the value of c is closer to 1, the attractor solutions of background become 
loose, toward the time when the vector field becomes heavy. More specifically, the equipartition of the vector energy density is not 
maintained to the sufficient level during the evolution. As is shown in FIG. [5] the full spectrum with Sip 5^ attains unacceptably large 
anisotropy for c > 1.25. 
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FIG. 4: Power spectrum P 7 , in the units of mt/Mp. The left panel shows the spectrum when c = 1.07, which achieves the 
observable level of statistical anisotropy, and the right when c = 1.5 in order to compare with FIG. [3] Results for some different 
values of £ are shown to visualize the directional dependence and its difference between the two cases. It is observed that the 
dependence on £ is flipped between the two panels: larger £ leads to larger P 1 (positive g,) on the left panel, and the opposite 
(negative <?») on the right. In fact, g„ monotonically decreases with increasing c, shown in FIG. [5] Here the mass of vector field 
is fixed to be mo = 1000 m v . 



Comparing the two panels, we observe that changing the value of the vector mass does not make a significant difference 
in the spectrum. The overall values of P 7 for the coupled case decrease by < 10% when mo is changed from 50 m v 
to 1000 Tji^p , but the angular dependence does not change: g* ~ —0.8 for both masses (as a comparison, ~ 0.15 
for the Sip = case). Although this value is undoubtedly excluded by observations, this comparison tells us that (i) 
the interaction between the vector and scalar through the vector kinetic and mass terms modulates the curvature 
perturbation significantly, that (ii) satisfying the conditions on the vector mass (M <§; a initially and M > a by the 
end of inflation) is not sufficient to attain the statistical isotropy, and that (iii) changing the vector mass does not 
affect the directional dependence of the power spectrum. 

The directional dependence of P 1 is sensitive to the value of c. FIG. [4] compares two cases of different values of c 
with the same vector mass mo = 1000 m^. For the purpose of comparison, the right panel takes the same choice of 
parameters as in the right panel of FIG. [3] (with two more sample values of £ and without Sip = case). In the left 
panel, the spectra for the case c = 1.07 are shown. For this particular value of c, the spectrum acquires the observed 
level of statistical anisotropy. It is interesting to see that the dependence on £ is opposite between these two cases. 
The spectral amplitude is larger for larger £ in the c = 1.07 case, and it is smaller for larger £ in c = 1.5: this is 
equivalent topositive for the former and negative for the latter. To quantify g*, we show its values as a function 
of c in FIG. [5] 10 It is found that g* does not change over the range of momenta. 

There are a few things to note. In every case, a nearly scale-invariant power spectrum is obtained, as exhibited in 
the figures. Secondly, acquires negative values in most of the parameter space, while it becomes positive for the 
small range 1 < c < 1.1. This implies that for c > 1.1, the larger pt is, the larger the spectrum, and it is peaked when 
p is completely orthogonal to the direction of the vector vev, i.e. £ — > 0. Thirdly, the interaction effect is maximal 
when £ is close to 0, as can be seen from FIG. [3] Since the momentum is completely orthogonal to the background 
vector vev at £ = 0, the larger amplitude for £ — > suggests that the interaction depends largely on px- Lastly, 
although the spectrum can achieve the observed anisotropy, the allowed window of c (g* is insensitive to mo) is very 
narrow, due to the strong dependence on c near g* ~ 0.3. This c dependence of the statistical feature of the spectrum 
is indeed real, resulting from the Sp — SA^ interaction. 



The angular decomposition of power spectrum, j7l| , is essentially a multipole expansion for small g*. When |g„| approaches to 1, the 
validity of this expansion starts to break down. The higher-order terms (higher than quadrupole) becomes non-negligible, and indeed, 
the decomposition ^71\ does not approximate the actual values well. We should treat with care the values of \g*\ ~ 1, but it is clear 
that such large anisotropy is excluded by the observations. 



17 




FIG. 5: The values of g, as a function of the coupling constant c, for mo = 1000 m v . For smaller values of Too, the curve would 
shift slightly upward, but the general feature would be unaffected. 



V. DISCUSSION AND CONCLUSIONS 



We studied the primordial inflation with a vector field non-minimally coupled to inflaton through vector kinetic 
and mass terms, and computed the power spectrum of curvature perturbations within the framework of curvaton 
mechanism. In this model pi), the mass term inevitably introduces a vector longitudinal mode, but it does not 
acquire a negative kinetic term and therefore is free from ghost instabilities. It is claimed in [TU] that a scale- 
invariant statistically-isotropic power spectrum can be obtained in this model. However, in their work the anisotropy 
of background evolution was neglected in the pure de Sitter universe, and the kinetic and mass terms were taken to 
be external functions of time and thus did not introduce any interaction effect due to the vector-inflaton coupling. In 
[39) . anisotropic background expansion was taken into account, and it was shown that the background attractor leads 
to the desired time-dependent functions, which we also showed in our language in Section [ii] However, the interaction 
effects were still assumed to be negligible. In this paper, we relaxed this assumption and did an extensive analysis in 
the linearized perturbation theory with the consistent inclusion of the vector-inflaton coupling. 

The non-trivial features arise due to the presence of the vector vev and its perturbations. The vev of vector field 
changes the dynamics of the system in three ways: (i) it modifies the background by slowing down the inflaton motion, 
(ii) it renders the direction of the vev a "privileged" direction, that expands differently from the other two, and (iii) 
it introduces interactions between the inflaton and the vector perturbations already at the linearized level. We have 
verified that the first two points do not make a dominating change. Namely, if we include (i) and (ii) (not (iii) yet), 
the anisotropy in the background expansion is small throughout inflation and is rapidly washed away after the vector 
field becomes heavy; also, the near scale invariance and statistical isotropy of the power spectrum can be attained 
(see dashed lines in FIG. [3]). The last point, however, generates significant effects on the spectrum. Although the 
scale invariance is not affected, the statistical isotropy is no longer a generic feature. We stress that this effect was 
disregarded in all the existing computation of the perturbations in this model. 

There are two parameters in the model that control the coupling ( 13 ) between the inflaton and vector fields: the 
inflaton-slowing parameter c, and the vector mass end value mo- While changing mo does not make any appreciable 
impact on the spectrum, the directional dependence of the spectrum is highly sensitive to the values of c. The measure 
of the statistical anisotropy, g„, is found to be a monotonically decreasing function of c and is mostly negative, except 
for a small range 1 < c < 1.1. This means that for c > 1.1, the value of spectrum increases as the direction of 
associated momentum approaches to the direction perpendicular to that of the vector vev. Also, as FIG. [3] suggests, 
the effect from SA^ — Sip coupling is strongest when the momentum is orthogonal to the vector vev. All these features 
can be justified by observing the action in the early- and late-time regimes. From (70), we see that only SA X mode 
contributes (at the linearized level) to the final value of the spectrum, so we look into the relevant part of the action. In 
the early time, we know M -C p, A ~ MA, and SA X ~ 5A y ~ Sip (order of magnitude, assuming small background 
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anisotropy) . Then the part of the action ( 29 1 which contains SA X simplifies to 

5 A y + h.cA -Pt \5A x \ 2 + °rVLVT (8A x 5A y + h. 



S. 



(2) early 
2dS 



dtd?ke 3a ' 



2 a 2 
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SA C , 



a p L p T 



b p 2 



f P 2 



A 



\Sip^ 5A X + h.cj 



(72) 



where . . . denotes the terms without 8A X . From (72), it is seen that initially the interaction term (second line) is 
negligible, as the initial vacuum is adiabatic, but in the time period when M ~ p, this term can make an O (1) 
contribution to 8A X , and thus to the spectrum. This term also illustrates that the interaction strength is proportional 
to c and pr, and also that it is unrelated to niQ. 11 In the late time, when the vector mass becomes M ^> p, all the 
directional dependence vanishes from the action, and so no anisotropy is further produced. Therefore, the statistical 
anisotropy in the spectrum is generated during M ~ p through the interaction between Sip and 5A X , and its strength 
is controlled by c and pr- This heuristic argument justifies the behavior of our numerical results. 12 

The vector field is a curvaton and needs to eventually decay to imprint a signature in the curvature perturbations. 
The curvature power spectrum is either |P 7 or |-P 7 , depending on whether p A or p T dominates the total energy 
density, but the features in the spectrum, such as scale invariance and statistical anisotropy, are directly inherited 
from P-y. While a nearly scale-invariant spectrum is a generic result of this model, statistical isotropy is not. We 
have shown that the observed value of the anisotropy parameter c/*, albeit controversial in itself, can be achieved in 
this model; however, the allowed parameter space is small (only near c = 1.07). This is in contrast with the results 
of [THIHO] where statistical isotropy was found to be a result of simply the classical evolution of / and m, and with 
the corresponding claim of [39 , in which the results for the perturbations of [T^l [2U] were used, implying statistical 
isotropy for any value of c. This discrepancy is due to the direct coupling between SA^ and Sip, which was neglected 
in [19, 20 , but should generically be present in any effective field theory realization of the mechanism. 

Throughout our analysis, we always assumed the initial equipartition of the vector energy density. In principle, 
we could compensate the unacceptable statistical anisotropy by adjusting the initial partition. However, since this 
partition is maintained during the evolution, it would still require fine-tuning. Thus it would not change our conclusion: 
the statistical isotropy of the curvature power spectrum is not in general attainable without finely-tuned parameters 
and initial conditions. 

Although this model is certainly not the only possibility to realize the vector curvaton scenario with varying kinetic 
and mass terms of the desired time dependence, we have studied a concrete, realistic example that consistently takes 
into account the vector-scalar interaction in the presence of the vev of vector field. Curvaton mechanism was first 
introduced in order to separate the generation of curvature perturbations from the detail of inflation dynamics and to 
relax the requirements for inflaton |18j . This mechanism is favorable in this sense, if the fields are minimally coupled. 
However, when multiple fields couple to each other in a non-minimal way, as required in |19U20) . the interaction non- 
trivially modifies the situation. In such models, calculation in the decoupled limit is not sufficient, and the consistent 
treatment of the whole system is crucial. 
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One might be concerned tha t mo should appear as there is A prese nt in the interaction term. However, the value of A (and A) is 
enforced by the attractor ^23^ . Since pa = ip 2 + D 2 ) and since |23| is independent of mo, we have \C\ = |D| <x — . Thus the 
interaction is not affected by the value of mo- 

As long as the attractor |23| is a good approximation, it is enforced that A <x \Jc— 1/c, which implies that the interation term in 
is oc \Jc — 1. Then it appears that as c — > 1, the effects from this term vanish and the statistical isotropy would be achieved. However, in 
the real evolution, the time variation of C and D (especially the latter) becomes appreciable as c — > 1, and the assumption of constant 
C and D is violated. Then the simple argument above does not apply; as c gets closer to 1, g* becomes positive and large. 
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Appendix A: Matrices in Early Time Evolution 



In this appendix, we show the complete expressions of the early-time matrices T, X and , appearing in (49 1 



As mentioned in a footnote in the main text, we did not antisymmetrize A for simpler expressions. However, 
this affects neither the initial conditions at the level of approximation we concern, nor the equations of motion, 
since antisymmetrization is simply to add total derivatives in the action. On the other hand, f2 2 is symmetric, by 
construction. 
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It can be seen from the above expressions that all the matrices are real and invariant under the parity p — > —p. 
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